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Distorted vortex lattice in a tetrahedral superconductor 
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Equilibrium shape and orientation of vortex lattice are studied for an s-wave tetrahedral superconductor 
in the vicinity of the upper critical field. The phase diagram, which includes transitions between rhombic and 
rectangular lattices, is constructed in the parameter space of the Ginzburg-Landau functional. The developed 
theory is applied to the heavy-fermion superconductor PrOs4Sbi2. In a wide range of parameters the shape 
of the vortex lattice is only weakly dependent on temperature. The neutron scattering measurements of the 
vortex lattice in PrOs4Sbi2 can be explained by a peculiarities of the tetrahedral symmetry group and are 
further supported by analysis of the appropriate band structure calculations. 



PACS: 74.20.De, 74.25.Qt 

Hexagonal vortex lattice predicted for ideal isotropic 
superconductors ^ |5j is perturbed in real materials 
by crystalline anisotropy. Anisotropic nonlocal cor- 
rections within the Ginzburg-Landau theory or in the 
London approximation are determined by the Fermi 
surface geometry (3j. Their effect may result in a 
hexagonal-to-square vortex lattice transition observed 
in the past in the superconducting borocarbides 0] 
[SI C]. Anisotropy of the Cooper pairs wave function 
'also gives rise to vortex lattice distortions. The line 
nodes of the dx2-y2 superconducting gap in the high- 
Tc cuprates favor, for example, a square vortex lattice 
|S1|H|. The flux line lattice (FLL) geometry provides, 
■thus, a combined insight into the Fermi surface and the 
brder parameter anisotropy. 

Filled skutterudite compound PrOs4Sbi2 has at- 
tracted significant attention in the past few years as 
a first example of Pr-based heavy fermion superconduc- 
tor with Tc = 1.85 K [iniiniinilTlini. At present, 
controversy remains about the symmetry and the struc- 
ture of the superconducting gap in PrOs4Sbi2. The 
low-temperature scanning tunneling spectroscopy ex- 
periments by Suderow et al. |15| clearly demonstrate 
]a finite superconducting gap over a large part of the 
Fermi surface. This finding is independently confirmed 
by an exponential low-temperature dependence of the 
nuclear relaxation rate 1/Ti l(i . though absence of 
the coherence peak may point to anisotropic pairing. 
Recently, small angle neutron scattering measurements 
have found a distorted vortex lattice in PrOs4Sbi2 at 
low fields and temperatures jJT]. The observed distor- 
tion was attributed by Huxley et al. to an anisotropic 
multicomponent superconducting order parameter with 
point nodes. PrOs4Sbi2 has, however, a rather unusual 



tetrahedral point group. This group contains three- 
fold rotations about cube diagonals but no four-fold ro- 
tations in contrast to the other cubic groups O and Oh- 
Investigation of the shape of FLL has not been done to 
our knowledge for such superconductors. In the present 
work we investigate the geometry of FLL in tetrahedral 
superconductors within the nonlocal Ginzburg-Landau 
theory for a conventional s-wave order parameter. 

The Ginzburg-Landau (GL) energy density func- 
tional for an s-wave superconducting order parameter 
^' derived from the BCS theory has the standard form: 



FGL = ^*i^){a + :Fu)^i^) + ^\^{^r+'^. (1) 

Here a = (1 — T/Tc), Tc is the transition temperature, 
and P = 7C(3)/(47rTc)2. The gradient terms 



•^n = ^ ^ H,, 

n—1 ii . . .i^n 



.H„ 



(2) 



are expanded into even powers of the operator H^ = 
—idi + (27r/$o)^ii 'i'o being the flux quantum. Expan- 
sion coefficients are expressed via the Fermi surface av- 
erages of the components of the Fermi velocity wp 'as 



K, 



^J5i(2-^)C(2n + l)(^^«..-F,)FS. (3) 



The tetrahedral point symmetry Th imposes certain re- 
lations between the gradient term coefficients: 



K^2 = K„ 



(4) 
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K,j.2y2 = Kx2y2 = Ky2^2^ K ^6 = Ky6 = K , 
Kxiy2 — Kyi;^2 — K^ix^, Kx2yi = K y2 ^4. = K ^2 . 
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Note, that for the tetrahedral group Kxiy2 ^ K^2yi. 

Let us now assume that an external magnetic field 
is applied along the l-axis. Considering only solutions, 
which are uniform along the field direction (n^^' = 0), 
the gradient terms are simplified to 



J'2 = K^ 



ui 

n 

n 



K,2y.{UlUl} , (5) 

K.,.y2{uiul} + K,2y.{ulnl}. 



In the above expressions { • • ■ } denotes a sum over all 
possible permutations of the gradient operators n,; and 

The upper critical field Hc2 is determined by the lin- 
earized GL equation, which can be conveniently writ- 
ten in terms of the Landau level operators a and a^: 
a — y^^oJ^7rH{Ilx — iHy). Rotation by angle if about 
5-axis transforms the Landau level operators into e~^'fa 
and e^'^a). After some algebra, the gradient terms ^ 
are expressed as 

T2 = h{2h + 1) , 

= [kio{2n^ -f 2n + 1) + ^44(0'' + a^^)] , 
= [A;6o".60 + fc62("-62a^ + a^^"-62) 

+k64:{n6Aa'^ + a^^fiQi) + keeia^ + a''^')] , (6) 

where h = 2-kHKj.2 /$o is a dimensionless magnetic field 
and 



kia=i{K,.+K,2y2) 12X^2 , 

fc44 — {K^i — 2>K,j.2y2)/2Kl2 , 

fceo = (2i^,6 + 3K,.y2 + 3K,2y.)/8Kl2 , (7) 

h2 = {Kx4y2 - K^2yi)/8K^2 , kee = -15^62 , 

fc64 = (2i^,6 - 5K,.y2 - 5K,2y4)/8Kl2 . 

The level-number operator h — d^a and its polynomials 
neo = {20h^ + 30n'^ + 40n+15), = (15n2-H63n + 45), 
^64 — (6". + 15) are invariant with respect to an arbi- 
trary rotation about the z-axis. The discrete rotations 
of the point crystal group are responsible for the ap- 
pearance of a" terms. In particular, the operators 
and break the four-fold rotational symmetry about z 
and discriminate between the x- and the y-axes. 

We use the standard procedure to determine the 
FLL geometry in the vicinity of the upper critical field 
PE]- Previously, such an approach has been applied 
for superconductors with tetragonal [HI 13 El E| j or- 
thorhombic [20], and hexagonal |2J| crystal structures. 
Solution of the linearized GL equations obtained from 
Eqs. IpJ and © is expanded in the Landau levels up to 
the sixth order: 




Fig.l: Geometry of the FLL unit cell as a function of C4 
and cg. Lines indicate corresponding values of angles, 
which parameterize a rhombic or a rectangular unit cells 

where /„ = (a''')"/o/Vn! and a/o = 0. In order to 
construct a periodic vortex structure at the zeroth Lan- 
dau level /o it is convenient to go from a laboratory 
frame determined by the crystal axes to a rotated frame 
such that a new i-axis points between a pair of nearest- 
neighbor vortices and is rotated by angle ip with respect 
to the crystal axis. The vector potential is chosen in the 
Landau gauge A = {—Hy, 0,0) and the periodic solu- 
tion with one flux quantum per unit cell is written as 

r I \ sr^ r • 2 < 27r . TriJ . 

/o(rj — y exp —mpm -) imx — ~ maa) 

(9) 

The basis vectors of the FLL are (a, 0, 0) and (ap, aa, 0) 
in the rotated frame, which satisfy the one-flux- 
quantum condition Ha? a — $0. 

The expansion coefficients in Eq. (jHJ are found for 
the eigensolution of the linearized GL equation as a per- 
turbation expansion in a small parameter a: 



V6, 45 

C4 = — A;44a , C2 = p 

4 2V2 



^62"^ , C6 = -\/5fc66a^ 



(10) 



* = AV , V - /o + C2e2*'^/2 +C4e^*'^/4 -l-cee^^'^/e , (8) 



Hip 



whereas the upper critical field is given by 

hc2 = -a - kioa'^ + (15fc6o - 2fc4o - 3kl^)a^ . (11) 

Neglecting the magnetic field contribution to the free 
energy /i^/Stt in the large-K limit we obtain for the en- 
ergy density 

{FGL)^X^{ria + Tn)^j)+X^^m') , (12) 

where (/) — (l/V) J dr^f{r). The quadratic term in 
the above equation is calculated as (|/oP)('i — ^c2)(l — 
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2k4Qa). Then, straightforward minimization of the en- 
ergy density ((T^ with respect to A yields 

^ 2/3(|^|4)/(|/o|2)2 ■ ^ ) 

The equihbrium distribution of the order parameter 
\E'(r) is found by minimizing the geometrical factor 

f3A = m')/{\fo\Y ■ (14) 

This generalized Abrikosov's parameter is a function of 
only three variables p, a, and (p. An explicit calculation 
yields 

Pa = v^^^ exp [27rip(TO^— n^) — 27Tcr{m'^+ n^)] Im.n, 

m,n 

(15) 

where summation goes over all integer and half-integer 
pairs {m,n). Function is defined by an integral 




where 



Piy) = 1 + C2e^'^H2iy) + c^e^'^H^iy) + cee'^'^Heiy) 

(17) 

and Hn{y) are the Hermite's polynomials. If we keep, 
for simplicity, only the terms, which are linear in c„, 
then the Abrikosov parameter is reduced to 

/3a « /3o + 4Re(c2e2*'^/32 + ae^'"^ (3i + cge^^'^/Se) , (18) 
where 

Pk = (/oVoA-)/(|/o|')' . (19) 

Function /3o (p, cr) is the standard energy parameter 
for an isotropic superconductor [Q, which is given by 
Eq. with Im]n = 1 p]- The other functions Pk are 
obtained from Eq. l(T3|l by substituting the correspond- 

mg Im'n- 

/W^ = i=(8.W-6W + ^) , (20) 
46)„ = -1= (le^V^n^ - 30. V^n^ + I wn^ - |) . 

Hexagonal vortex lattice with p = 1/2 and a = V3/2 
corresponds to the absolute minima of functions Po (p, a) 
and PeiPjCr). The other two functions Pi and P2 tend 
to stabilize a square and a distorted triangular lattices. 




Fig.2: Expansion of the right part of fig.l. The shaded 
region is composed by parabolas ce{T) oc C4{T)^; the 
black line is defined by cg — 25c4 

respectively. Their competition yields a rich phase di- 
agram of a tetrahedral superconductor. In finding the 
lowest energy vortex lattice we have kept also additional 
terms in Pa, which are proportional to high-orders of c„. 

Equilibrium form and orientation of the vortex lat- 
tice in a tetrahedral superconductor should be obtained 
by minimizing /^^(p, cr, (/?) for different sets of the pa- 
rameters a and kij in the GL functional. The task is 
considerably simplified if the expansion of tp in Eq. (|HJ| is 
rapidly converging, e.g., near Tc or for a superconductor 
with a weak anisotropy. In such a case, a set of the ef- 
fective parameters can be reduced to C4, which is deter- 
mined by the lowest order anisotropy in quartic terms, 
and cg, which quantifies the x—y discrimination intro- 
duced by the tetrahedral symmetry. Since C2 oc cg in the 
leading order in a, we fix C2 = — 3/(2vT0)c6 « — 0.47c6. 
The energy parameter Pa is, then, numerically mini- 
mized with the conjugate gradient method. The pa- 
rameter space is further restricted to Cg > since a 
change of the sign of cg corresponds to a rotation of the 
coordinate frame hy if = 90°. 

Our main results are presented in fig^ where the 
geometry adopted by the FLL is plotted for different 
values of C4 and cg. The considered region in the pa- 
rameter space is divided into two parts corresponding to 
highly symmetric vortex lattices separated by a transi- 
tion region colored in black. In the main area, the whole 
(c4,cg) plane apart from the lower left corner, the unit 
cell has a rhombic shape with the longest diagonal par- 
allel to the y-axis. The apex angle of a rhombus varies 
from 50° to 90° as indicated by thin dotted lines. 
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Since C4 ~ 0{a) and cg = 0{a^), an evolution along 
the Hc2{T)-cuTye follows approximately a parabolic 
path in the (04, C6)-phase diagram, which starts at the 
origin C4 = and cg = 0. For a certain range of pa- 
rameters, the apex angle along such an Hc2{T)-path 
does not change significantly with temperature. For ex- 
ample, 9 remains between 70° and 80° in the shaded 
area of figEl composed by a set of parabolas given by 
Ce/cl = — 16V5fc66/6fc|4 with 

-5 < kee/kl^ < -0.6 . (21) 

The reason for such a weak dependence is a cancellation 
of two opposite tendencies determined by (32 {p, c) and 
/34(yO, a), see Eq. ((T^ . If the crystal would have the cubic 
or the tetragonal symmetry, cq is zero and the apex an- 
gle increases almost linearly between 60° and 90°, with 
a phase transition into a square lattice afterwards. The 
present investigation shows that the tetrahedral symme- 
try introduces a fundamental difference via a sixth-order 
gradient terms. 

In the lower left corner of the phase diagram, for 
C4 < and small ce, the shape of the vortex unit cell 
changes via a second-order transition from a rhombus 
with the longest diagonal oriented by 45° from the x- 
axis to a rectangle with the longest side along the y-axis. 
When C4 decreases, the small apex angle of the rhom- 
bus goes from 60° to 90°. The two regions on the phase 
diagram are separated by a transition region (shown in 
black), where the unit cell of FLL has no mirror sym- 
metry. FigO illustrates that a transformation of the 
unit cell corresponds to a shear instability of the vor- 
tex lattice: vortex chains in the rectangular lattice slide 
along one of the directions such that the lattice turns 
into a centered rectangular (rhombic) lattice with the 
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Fig. 3; Expansion of the left side of fig.l: transition re- 
gion between the vortex unit cells of rhombic and rect- 
angular shapes 



Table L Fermi surface averages of different com- 
binations of the components of the Fermi veloc- 
ity {vp^VpyVp^)FS for the 49th band in units of 
(10^ m/s)''-+'+'". 



klm 


a-sheet 


7-sheet 


combined 


200 


2.19 


1.46 


1.61 


400 


14.7 


6.36 


8.11 


220 


1.51 


1.63 


1.60 


600 


127. 


45.0 


62.2 


420 


5.69 


3.74 


4.15 


240 


4.49 


8.75 


7.85 


222 


1.03 


0.63 


0.714 



same volume of the unit cell and the same side ratio. 
This transformation goes via two second-order transi- 
tions for C4 < -0.06, whereas for -0.055 < C4 < -0.048 
the transition is of the first order. As a last remark, we 
would like to point at a peculiar possibility for a tetra- 
hedral superconductor with negative C4. While moving 
along the Hc2{T) line towards low temperatures, the 
vortex lattice first changes from a triangular one into 
a rectangular lattice and then transforms back into a 
distorted triangular lattice. 

Finally, we shall apply the above results to 
PrOs4Sbi2. Topology of the Fermi surface of this 
material has been studied by the de Haas-van Alphen 
measurements and compared to the results of the 
LDA-j-C/ band structure calculations by Sugawara et 
al. |22[ 05| . The Fermi surface is composed of three 
sheets: one of the 48th band and two of the 49th band. 
The contribution to the density of states from the 48th 
band is relatively small: iV48(0) ^ 0.04iV49(0). Besides, 
only a small Fermi velocity in the 49th band is capable 
to account for a large value of the upper critical field 
at zero temperature i?c2(0) = 2.2 T ^Hl- Accordingly, 
we assume that the active band, which is responsible 
for superconductivity in PrOs4Sbi2, is the 49th band, 
whereas the 48th band plays a passive role and may 
have a smaller gap amplitude as suggested by Seyfarth 
et al. [m. 

The band structure data for 195 k-points in the ir- 
reducible part of the Brillouin zone of PrOs4Sbi2 have 
been previously used in the comparison between LDA 
results and de Haas- van Alphen data [221 ■ We inter- 
polated between these data points using appropriate 
lattice harmonics and, then, calculated numerically the 
Fermi surface averages for the 49th band using a much 
smaller mesh in the momentum space. The combined 
averages have been obtained by a sum of the two contri- 
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butions weighted according to partial densities of states 
N^g/Njg ~ 0.27. The results are summarized in Table 
I. Using Eqs. Q and ([TJ we find for the dimensionless 
gradient constants fc44 = —0.29, fcge — —0.35 such that 
^66/^44 = —4.3 and 

C6(T)/c4(T)2 « 25.0 . (22) 

The above relation is illustrated in figElby a solid line. 
For C4(T) > 0.02, which corresponds to \a\ > 0.12 and 
T < 0.9Tc, the vortex lattice has a rhombic unit cell 
with a nearly T- independent apex angle = 70°. For 
a rhombic lattice this angle coincides exactly with the 
angle between two reciprocal lattice vectors. The lat- 
ter angle has been measured in the neutron scattering 
experiment ^2] and is equal to 6* « 75 ± 5°. Thus, 
there is a good agreement between our calculation and 
the experimental data on PrOs4Sbi2. We emphasize 
again, that in a tetrahedral superconductor such a sta- 
ble distortion of a hexagonal vortex lattice arises due 
to a competition between several anisotropic gradient 
terms. When comparing the above results with ex- 
perimental data one should, of course, bear in mind 
that applicability regions are somewhat different. For- 
mally, our calculation has been restricted to the vicinity 
of Hc2, but is valid in a wider range of applied fields 
H > 0.2Hc2 in high-K materials such as PrOs4Sbi2 
(k Ri 29). The temperature range is in the GL regime 
O.GTc < T < 0.9Tc, whereas the neutron measurements 
have been performed for T < 0.45Tc. Further investi- 
gations of the shape of the vortex lattice in tetrahedral 
superconductors in the London limit at low tempera- 
tures would be, therefore, useful. 

A different confirmation of our analysis can be pro- 
vided by angle dependence of the upper critical field 
since the amplitude of the modulations are related to the 
anisotropy of the Fermi surface El ESI- ^'^^ example, 
a perturbation calculation yields the main contribution 
to the [001]-plane four-fold modulation: 

^ = 1 + ^-^44 (cos 40-1), (23) 

where (j) is angle between an applied magnetic field and 
the i-axis. 

We are indebted to V. P. Mineev for stimulating dis- 
cussion and to H. Harima for providing us the band 
structure data for PrOs4Sbi2. 
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